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Abstract. We investigate correlations of current at contacts and voltage 
fluctuations at voltage probes coupled to interferometers. The results are 
compared with correlations of current and occupation number fluctuations at 
dephasing probes. We use a quantum Langevin approach for the average 
quantities and their fluctuations. For higher order correlations we develop a 
stochastic path integral approach and find the generating functions of voltage 
or occupation number fluctuations. We also derive a generating function for 
the joint distribution of voltage or occupation number at the probe and current 
fluctuations at a terminal of a conductor. For energy independent scattering we 
found earlier that the generating function of current cumulants in interferometers 
with a one-channel dephasing or voltage probe are identical. Nevertheless, the 
distribution function for voltage and the distribution function for occupation 
number fluctuations differ, the latter being broader than that of former in all 
examples considered here. 



PACS numbers: 73.23.-b, 03.65.Yz, 72.70.+m 
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1. Introduction. 

The quantization of charge and the diffraction of quantum waves lead to fluctuations 
in the transport state of electrical conductors even in the limit of absolute zero 
temperature pQ . These fluctuations provide additional information not available from 
conductance measurements alone. This has led to considerable theoretical O [3j SI E] 
and experimental efforts [6] [71 02 [10] to characterize fluctuations not only on the 
level of noise but including all higher order cumulants. These efforts have been focused 
almost exclusively on the statistics of transferred charge [UOH], i-e- on the cumulants 
of currents measured at a contact of a conductor. Clearly, there are however also 
other quantities of interest. For example, fluctuations of the charge inside a conductor 
[Til [T2l [13] (as opposed to transferred charge) are closely connected to electrostatic 
potential fluctuations in a conductor [14] and are thus essential in determining the 
long range Coulomb interaction of two nearby conductors. Still another example 
is fluctuations of voltage. The statistics of voltage have been investigated for a 
mesoscopic conductor in series with a classical impedance [TS], Qj)] and for networks of 
two-terminal conductors [17] , 

In this work we consider voltage probes which connect to a conductor within 
its phase-coherent volume [18]. We determine the fluctuations of the voltage at the 
probe [UlEOlET] and determine the joint probability distribution between the voltage 
fluctuations and the charge transferred into a contact of the conductor. Voltage probes 
are real elements of mesoscopic conductors and are used to gain information of the 
interior state of the system. We compare the fluctuations at a voltage probe with 
the fluctuations at a dephasing probe. This latter probe conserves not only charge 
but also energy [22] [23] . At a dephasing probe the occupation number fluctuates [2_0j . 
Voltage and dephasing probes provide a simple means to introduce incoherent events 
in an otherwise quantum coherent conductor [181 123j . 

We examine these correlations with the help of a (quantum) Langevin approach 
|24[ 125] . In a second step we derive a generating function for the correlations using a 
semi-classical path integral approach [26] [27] . On the level of conductance and noise 
the Langevin approach and the saddle point solution of the path integral approach 
are shown to agree. We have earlier demonstrated that for the statistics of transferred 
charge both a single channel voltage probe and the single channel dephasing probe 
lead to the same counting functional if the scattering matrix of the conductor is energy 
independent in the range of applied voltages [26] [27] . Interestingly, we find that this 
equivalence does not hold for the joint distribution of voltage (or occupation number) 
fluctuations and transferred charge. 

The characterization of the transport fluctuations not only in terms of transferred 
charges but also of internal properties is a central motivation of this work. Voltage 
fluctuations at a probe test carrier fluctuations that can still be evaluated in an 
approach that considers samples with terminals only. It is particularly interesting for 
conductors which exhibit quantum interference since voltage and dephasing probes 
are phase breaking elements. A carrier that enters a probe is eventually replaced by a 
carrier that is re-emitted into the conductor with a phase (and in the case of a voltage 
probe, an energy) that is unrelated to the exiting carrier. This carrier is detected 
in the voltage probe and thus does not contribute to certain interference processes. 
Phase breaking and which path detection are generally closely linked phenomena [28J . 
This provides another motivation to investigate the fluctuations at the probes (the 
which path detectors) and their correlations with the transmitted current. 
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2. Correlations: Langevin approach. 

The Langevin approach is a convenient tool to investigate current and voltage 
correlations in mesoscopic conductors [24] [25]. This has been demonstrated in 
numerous papers, see e.g. Refs. [TJ [TH] 1301 1311 132] • Importantly, the Langevin 

approach also allows one to obtain information on the fluctuation properties of a probe 
connected to the mesoscopic conductor. In this section we investigate the average 
current and potentials as well as the auto- and cross-correlations. It is possible to 
extend the Langevin approach successively to higher order correlations, giving rise 
to cascade corrections [5]. In section [3] we instead directly determine the generating 
function of the correlations, making use of a stochastic path integral approach 0]. 

We consider a multi terminal quantum coherent conductor first without being 
connected to a probe. The current through the conductor fluctuates. When 
the conductor is embedded in a zero- impedance external circuit, all contacts are 
characterized by Fermi distribution functions which are constant in time. Only bare 
current fluctuations take place due to scattering in the conductor, described by a 
unitary scattering matrix S. The total current in contact a can be expressed as a sum 
of the average current and the current fluctuations [24J [25] 

Ia = (I a )+SI a = ^G afi Vff+SI a , (1) 



where Vg is a voltage applied to reservoir fj, and G a Q the conductance matrix 

G Q/3 = ^ J dE [N a S aB - T a0 ]. (2) 

The equilibrium Fermi function is called /, the number of channels in lead a is N a , 
and T a g — T a g(E) — Tr[Sl l g(E)S a g(E)] denotes the total probability that a particle 
is transmitted from terminal (3 to terminal a where the trace runs over the transport 
channels. 

Like the conductances, the noise correlations can be expressed in terms of 
scattering matrices. The zero frequency noise and correlations are C a g = (SI a SIg). 
Note that the definition of the noise as given for example in Ref. 1 usually differs by a 
factor of 2. The noise can be expressed as C a g = C^L + C^g, a sum of an equilibrium 
contribution CJg stemming from noisy incident beams in one contact proportional to 
/q(1 — fa) j and a transport contribution C£a which depends on Fermi functions of 
two different reservoirs / 7 (1 — fs) (with 7 7^ 5). Following Ref. [25] we find 

C Z = \ J dE(26 a0 N a f a (l - f a ) - 

Tp a f a 0—f a )-T aP f p (l-fp)), (3) 

- fs) + fsO- - A) - 2/ Q (l - /„)) . (4) 

In the case when temperature is negligible compared to the applied bias, the 
contribution C^L vanishes and the term C^g represents the pure shot noise. 
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2.1. Voltage probe. 

A voltage probe is a large metallic contact connected to the conductor. The potential 
V p at the contact is left floating and there is no current drawn at the probe. The 
potential V p = V p (t) exhibits fluctuations on the time scale Td, which is given by the 
i?C-time of a classical circuit. The fluctuations originate from the response of the 
potential to the injected charges: Incoming charges make the potential V p rise. This 
leads in turn to an increase in the outgoing current which reduces the potential again 
and so on. This mechanism is formally expressed by vanishing current and current 
fluctuations at the probe for frequencies smaller than X/tj,: 

(I P ) = 0, AI P = 0. (5) 

The voltage probe is described by an equilibrium Fermi function f P (V p ), where 
thermalization by inelastic scattering is assumed to be much faster than the delay time 
Td, and where the temperature is determined by the lattice temperature. Therefore, 
a particle scattering via the voltage probe looses not only its phase memory but also 
changes its energy. 

The potential fluctuations at the probe give rise to additional fluctuations of the 
current in the terminals [El [20] . The potential at the probe can be separated into a 
constant and a fluctuating part, V p = V p + AV p . Following Eq. JT]) the average current 
(I a ) and the total fluctuations AI a are in linear response given by 

(I a ) = G ^pVp + G ap V pi (6) 
/Mp 

AI a = G ap AV p + SI a . (7) 

Both the constant value V p and the fluctuating part AV P of the potential at the probe 
are determined by the condition Eq. ([5]) . Throughout this paper we consider that only 
one of the contacts (taken to be a = 1) has an elevated potential eV and all other 
current terminals are grounded. Solving Eqs. ^ and (J7|) with a = p, the potential at 
the probe is found to be V p — ~G p iV/G pp , and AV P = —SI p /G pp . Using this we find 
for the auto-correlations of the potential and the cross correlations of the potential 
with the current in terminal a 

(AV p 2 ) =§F> (8) 

(AV p AI a ) = (811) - (SI p SI a )) IG PP = QsggL _ <k. (9 ) 

\ Wp / ^pp ^pp 

These expressions can be evaluated for a given scattering matrix using Eqs. @, §5§ 
and ([3]), where only the constant Fermi functions enter, i.e. f p — f p (V p ). For higher 
order correlations the fluctuations of the Fermi functions have to be taken into account. 
As pointed out above, this leads to noise of noise and results in a cascade corrections 
similar to the one known from the quasi-classical fluctuating Boltzmann equation [5] . 

2.2. Dephasing probe. 

Dephasing probes are conceptual tools to describe quasi-elastic dephasing [551 • In 
the case of a dephasing probe, a particle entering the probe is incoherently re-emitted 
within the same energy interval [E,E + dE\. It looses its phase, but the change 
in energy dE is much smaller than the external bias or the temperature. Because 
scattering in each energy interval is independent, the distribution function n p (E) is 
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-unlike the Fermi function describing the voltage probe- a highly non-equilibrium 
distribution function. It shows independent fluctuations in each energy interval, 
n p (E) — n p (E,t), which occur on the time scale tj, the delay time of the probe. 
In a dephasing probe, the current per energy interval i p (E) and its fluctuations are 
conserved up to the frequency 

(i p ) = 0, Ai p = 0. (10) 

Similar to the potential of the voltage probe, the non-equilibrium distribution function 
is written as a constant and a fluctuating part f p = n p + An p . All other terminals are 
still characterized by constant equilibrium Fermi functions. The average current per 
energy (i a (E)) — (i a ) and the total current fluctuations per energy Ai a {E) = Ai a 
are given by 

(i a ) = ^gapfp + g ap n p , (11) 

Ai a = g ap An p + Si a . (12) 

In analogy to the expressions for the energy integrated current I a = j dEi a (E), an 
energy dependent conductance matrix g a p = g a p(E) can be defined, 

g a p = j\N a 8 a [3 - T a p]. (13) 

Also for the zero-frequency noise we find an energy dependent function from C a p = 
j dEc a p(E). Again the function c a p = c a p(E) — (5i a 5ip) splits into an equilibrium 
and a transport contribution 

^=c e : P + c %, (u) 

where c^g and c^g are the integrands of Eqs. (HJ) and ((4|) respectively. c a q p contains 
the noise due to a noisy incident beam in one contact, while cvg involves particles 
from two different contacts. 

Using Eq. (fT0|) the noise and correlations at the probe are [2Ql [33] 

(An*) =(^)/4 = ^, (15) 

9ap /r-2\ /c- c- \ \ / _ 9apCpp Cp a 



(An p Ai a ) = ^ (Si*) - (5 lp 8 la ) /g pp = - ^. (16) 

V 9pp / 9 PP 9pp 

Here we can see the first consequence of the different electron distributions (see Fig. [TJ 
in the voltage and dephasing probe. Although Eqs. ((SJ) and (fl~5|) have exactly the 



ey, eV 



b) 



n p + 



eV 



Figure 1. a) The Fermi distribution function of the voltage probe takes on only 
the values zero or one (at zero temperature). The position of the step is given by 
the potential at the probe, b) The distribution function of the dephasing probe is 
a two-step function and has a value between zero and one in the interval [OjeV]. 
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same structure, the evaluation leads to different results. The reason is that in the case 
of a voltage probe the contribution C p ^ to the noise vanishes at zero temperature, 
while for a dephasing probe, the corresponding expression is non-zero even at zero 
temperature since the factor n p (l — n p ) is finite in the energy interval [0, eV] of interest. 
Interestingly, in Eqs. ([9]) and (fT6)l . inserting the expressions for the conductance and 
correlators from Eqs. ©-Q, it follows that the terms proportional to n p (l — n p ) 
drop out. Consequently, for energy independent scattering, the energy integrated 
correlations (AV p AI a ) /V equal the correlations (An p Ai a ) . 

2.3. Correlations and interference. 

It is interesting to investigate how and if interference appears in the fluctuations of 
the potential and the occupation number of the probe or in the correlations with the 
current in one of the leads. To this end we evaluated Eqs. (jHJ), {2J and (|15|. (|16p for 
four different setups: a beam splitter and three interfering systems, the Mach-Zehnder 
interferometer, the double barrier and the triple barrier. 

These examples also illustrate an interesting feature of the different nature of 
voltage and dephasing probes. It was shown in Refs. [26, 27J that the auto- and cross 
correlations (actually all cumulants) of the current in a conductor coupled to a single 
mode probe are independent of whether the probe is a voltage or a dephasing probe for 
the case that scattering is energy independent. Interestingly, this does not hold for the 
fluctuations of the probe: there is a clear difference between the potential fluctuations 
of a voltage probe and the occupation number fluctuations of the dephasing probe. 

All the examples are quantum Hall bars of different complexity subject to a high 
magnetic field. Current is transported by an edge state at filling factor one, where 
carriers move in one specific direction along the edge of the structure. A quantum point 
contact in such a structure forms a beam splitter: a particle is either transmitted and 
continues along the same edge, or it is reflected and moves into the opposite direction 
along the opposite edge. The beam splitter is the smallest unit that the interfering 
structures are composed of, we chose it here as a simple and instructive example 
to show how a voltage or dephasing probe behaves. For simplicity we consider the 
scattering matrix of all structures to be energy independent. In the following we 

P 



A V 




Figure 2. A quantum point contact in the quantum Hall regime forms a beam 
splitter with transmission probability T. A perfectly coupled probe is attached 
along the upper edge. 
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Figure 3. A Mach-Zehnder interferometer with its upper arm coupled to a 
dephasing or voltage probe with coupling strength e. An Aharonov-Bohm flux 
generates a phase 



present briefly the setups of the four examples. The results are collected in table Q] 
and table [H 

• The beam splitter [Ml [35] is formed by a quantum point contact with 
transmission probability T, situated between two terminals 1 and 2, as shown 
in Fig. [2] A probe -either a voltage or a dephasing probe- is perfectly coupled to 
the setup, such that every particle transmitted through the beam splitter enters 
the probe and moves then from the probe into contact 2. 

• An electronic Mach-Zehnder interferometer (MZI) is shown in Fig. [3] It 
consists of two arms connected to four electronic reservoirs 1 to 4 via beam 
splitters A and B. The transmission (reflection) probabilities of the beam splitters 
are Ta and Tb {Ra and Rb) respectively. Interference occurs, because the 
electrons have two alternative paths to propagate through the interferometer 
between beam splitter A and B. An Aharonov-Bohm flux threads the two arms, 
and the different vector potentials in the two arms lead to a phase difference 
$. This difference creates a characteristic flux-periodicity in the interference 
pattern, the Aharonov-Bohm effect. A voltage or dephasing probe is attached to 
the upper arm. Particles entering the probe loose their phase and, therefore, the 
interference is reduced, controlled by the coupling strength e. Thus the Mach- 
Zehnder interferometer is useful as a conceptually simple interferometer and has 
been used in theoretical discussions of dephasing [36 , 37, 38, 39] 00]. Such a setup 
was recently realized experimentally [HI S3] and has since generated further 
works [HI 25J US] H7] . If the length of the two interferometer arms is not equal, 
the transmission through the interferometer will be energy dependent due to the 
geometric phase a particle acquires. For finite temperature, interference effects 
will be averaged out. In order to determine purely the effect of dephasing and 
voltage probes, we concentrate here on the ideal case of an equal arm lengths 
interferometer. 

• The double barrier, shown in Fig. 21 is a two terminal conductor with two 
beam splitters @8l|49]. In the region between the barriers (the dot) a particle can 
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perform multiple loops and undergo resonant scattering. For each loop in the dot, 
the electron picks up a phase We consider zero temperature and an energy 
independent phase. The interference is here not a consequence of two spatially 
distinct paths but of multiple possible paths through the closed orbit inside the 
dot. A probe is coupled to the dot with coupling strength e. Since particles 
that scatter via the probe loose their phase, coherence is partially destroyed. For 
perfect coupling, e = 1, all particles entering the dot pass on into the probe, then 
the system is incoherent and in the sequential tunneling regime. 

• In the triple barrier setup, an additional barrier in series is inserted into the 
double barrier setup, see Fig. [5l Now there are two dots (EH [5TJ [52] with two 
different phases $l and Qr. Resonances can occur in either one of the dots or 
through both dots simultaneously, rendering the scattering matrix considerably 
more complex. The limiting case i?2 — *■ however represents the double barrier 
with + <f>^ = We consider here a probe coupled to the left dot, thus only 
the interference in the left dot is destroyed by the probe. 

In table Q] the results for the auto-correlations of the voltage (AVp) and of the 
occupation number (An^) are listed. The different electron distributions of voltage 
and dephasing probes, as pointed out already in section \2l2\ become apparent: the 
fluctuations of the occupation number n p (E) in a dephasing probe are larger than the 
potential fluctuations of the voltage probe. This is most easily seen in the example of 
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Figure 4. The double barrier quantum Hall interferometer with a resonant state. 
A probe couples to the resonant state with coupling strength s. $d is the total 
phase acquired in one loop. 



P 




Figure 5. The triple barrier with two resonant states. A single probe is coupled 
to the left resonant state. 
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Table 1. The mean squared fluctuations of the potential and the occupation 
number for the four examples. The transmission probabilities of all point contacts 
in the double and triple barrier structures are equal, T = l—R. Interference effects 
are present only in the latter two. 





(AV 2 ) -2- 

/ hV 




Beam splitter 


T(l-T) 


2T(1 - T) 


Mach-Zchnder 
interferometer 


*=ZT A (1-T A ) 


ZT A (1-T A ) 


Double barrier 


R{{l+R 2 )(2-s)+AR^/T : ^e cos <S> d ) 


2R(1 + R 2 {l-e)+2R^T^I cos $ d ) 




e(l-R){l+R) 3 


Triple barrier 


f^( J R,$ fl )[(l + 3i? 2 )(2-e) + 

4i?(l + i?)Vl-ecos$L + 
4i? 2 Vl - ecos($ L - $ H ) + 
2i?(l + i?)(2-e)cos$ fl + 
4i?Vl - ecos($ L + 

with 

-p(r> if.„\ _ l + R 2 +2Rcos<S> R 


?fT{R^ R )[l+m 2 -2R 2 e+ 

2i?(l + i?)v / l-ecos$ L + 
2i?Vl-ecos($ L -$ 7? ) + 
2i?(l + i?(l-e))cos$ jR , + 

2Ry/l -£C0S($ L + 


j-{it,^ R ) - {1+R+2Rcosi , R yi 



the beam splitter, where the fluctuations of n p (E) are twice as large as the fluctuations 
of V p . Note that occupation numbers at different energies are uncorrelated but that 
the fluctuation spectrum, in the limit of zero temperature, and in the limit that the 
energy dependence of the scattering matrix can be neglected, is independent of energy 
for < E < eV. 

For finite coupling, the fluctuations are proportional to f je and diverge for e — * 
for both voltage and dephasing probe. This means the fewer particles enter the probe, 
the stronger are the fluctuations. Note that the differences in the fluctuations of 
voltage and dephasing probes are of order one compared to 1/s and, therefore, become 
less important for small e. (For the beam splitter, the table includes only the strong 
coupling result e = 1. The above statement applies also to the beam splitter with a 
weakly coupled probe). 

The noise of the probe in the MZI is independent of the phase $. This is 
not surprising since particles enter the probe before they could interfere. This is 
in accordance with Ref. [36] which demonstrates that charge density fluctuations in 
an MZI arm are independent of phase. In the double barrier, particles enter the probe 
after performing multiple loops inside the dot and thus the fluctuations depend on the 
phase $d picked up. In the triple barrier different paths containing different resonant 
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Table 2. Current voltage correlations {AV p AI a ) and current occupation number 
fluctuations {An p Ai a ). For energy independent scattering these two correlations 
are identical. 





(AV p AI a ) ± = (An p At a ) ± 


Beam splitter 


±T(1-T) 
for a = 1 (+) and a = 2 (-) 


Mach-Zehnder 
interferometer 


T A (1-T A )(1-2T B ) fora = 3 


Double barrier 


tW for a = 2 


Triple barrier 


2i? cos $ b] for a = 2 



passages are possible and the noise contains a combination of $_r, $l + $_r and 

Table [2] shows the correlations of the voltage fluctuation measured at a probe with 
the current transferred into a contact. As usual p denotes the probe and a stands for 
any of the contacts of the conductor. As pointed out at the end of subsection 12.21 the 
energy integrated correlations (AV p AI a ) /V equal the correlations (An p Ai a ) which 
characterize an energy interval in the window opened by the transport voltage eV. 

The example of the beam splitter shows that the correlation can have either sign, 
depending on whether we consider the correlation with the transmitted (a — 2) or the 
reflected (a = 1) particle current. The dependence on the transmission probability 
directly reflects the shot noise generated at the quantum point contact. 

The correlations in the MZI are independent of the phase <&. The correlations are a 
signature of incoherent processes: only carriers that enter the probe give rise to voltage 
fluctuations and only carriers that leave from the probe to the contact are correlated 
with the voltage at the probe. Note, that the additional beam splitter B the particles 
have to pass after leaving the probe renders the correlations cubic in the transmission 
of the beam splitters (in contrast to the beam splitter setup of the first example). 
Curiously, for Ta = Tb = T they coincide with the third current cumulants [H [33] 
of a beam splitter with transmission T. In contrast to the simple beam splitter (first 
example) here the correlation changes sign as function of the transmission probability 
of the second beam splitter Tb- For Tb = 1 we have the same situation as for the 
beam splitter for a = 2, and Tb = corresponds to the beam splitter result for a = 1. 

Surprisingly, also the correlations for the double barrier are independent of the 
phase The fact that no oscillating contribution appears in the correlations is 
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however a feature of this particular chiral setup. For non-chiral Aharonov-Bohm 
rings, for example, the current- voltage correlations do depend on the phase 4^ [55]. 

In the triple barrier a probe is only connected to the left dot. The correlations 
of the triple barrier contain the phase due to interference in the right dot but are 
independent of $l . This is consistent with the fact that the correlation for the double 
barrier is independent of phase. 

All correlations shown in table [2] are independent of the coupling strength 
of the probe and moreover are independent of the phase directly adjacent to the 
probe. Further research is needed to characterize the examples where current-voltage 
correlations are phase- and coupling-strength independent and the examples where 
this is not the case [551. 



3. Full counting statistics of the probe: the stochastic path integral. 

An elegant description of the zero frequency transport is the full counting statistics 
(FCS) 2, 3, 5 J which permits to obtain not only the average current and noise but the 
full distribution of charges transmitted through a conductor during a measurement 
time r. For an M-terminal conductor without a probe the distribution function is 
denoted by -P(Q), where the vector quantity Q = (Qi, Q%, ■ ■ . , Qm) describes the 
charge transfered into each of the M terminals. P (Q) can be expressed in terms of 
the cumulant generating function 5(A), where A = (Ai, A2, . . . Am) are the conjugate 
variables to Q: 

P(Q) = / dAe s ^~ lA - Q , (17) 



5(A) = ln^P(Q)e ,;A - Q . (18) 
Q 

The sum and integrals run over all elements of the vector,/ dA = (2tt)~ m J d\i . . . d\M 
and = 5Zq 1 q m ■ Probability conservation leads to the normalization of the gen- 
crating function 5(0) = 0. 

For a long measurement time t, the transmitted charge into a contact a is 
proportional to r: Q a ~ rl a . Then the zero frequency current cumulants are obtained 
by taking derivatives of the cumulant generating function with respect to the counting 
variables and evaluated at A = 0. The average current (I a ) and the auto- and cross- 
correlations C a f3 are given by 

, x e dS e 2 d 2 S 

(la) = "T "TTT j Cap = — a > a, ■ (19) 

When a probe is connected to the conductor, particles can enter and leave the 
probe. However, for every particle entering the probe, a particle is re-emitted after 
a delay time Td as described in sections [2~T1 and l2~2l The functions P(Q) and 5(A) 
have then to be determined under the condition that no charge is accumulated in the 
probe. The transport statistics under this constraint was developed in Refs. [27] 
using a stochastic path integral approach. The findings are that the model of voltage 
and dephasing probes are perfectly equivalent for the case of one single-channel probe. 
For multi-channel or multiple probes the transport statistics of voltage and dephasing 
probes differs. 

Here we address the question of whether information about the statistics of 
the probe itself can be obtained. Usually the quantities investigated with the FCS 
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Figure 6. a) The current as a function of time at a terminal. The charge 
transferer! into a terminal is absorbed immediately. For a one-channel contact 
only one charge per time step r B = h/eV can enter, b) Charge on a voltage 
(dephasing) probe as a function of time. The charge on the probe accumulates 
and fluctuates on the time scale t&, the delay time of the probe. The figure is for 
Td > Ts ■ 



arc the number of charges Q = J Q r dtl(t) counted during a measurement time r. 
Consequently, we now need to consider the statistics of the quantities J dtV p (t) and 
J dtn p (t). More precisely we are interested in for example A j dtdt' (V p (t)V p (t')} 
or - J dtdt' (n p (t)n p (t')) . For the voltage probe we can define a normalized phase 
4> p proportional to the time integrated voltage and for the dephasing probe a time 
averaged occupation number n p = n p (E), 

^ = ^11 dtv ^ = ihL dtVp(t) ' (20) 

n p = -[ dtn p (E,t), (21) 

T JO 

where N — eVr/h. These are the accessible quantities describing the probe, they 
are not numbers of particles but phases. This is a consequence of the fact that 
charges transfered into the terminals are absorbed, but the charge on the probe is 
conserved. Fig. [5] illustrates this difference. Both <p p and n p vary from measurement 
to measurement, because the injection of particles into the probe that determines the 
time averaged potential and the time averaged occupation number is a probabilistic 
process. The phase 4> p is proportional to the total charge on the probe integrated 
during the measurement time r: <p p = c e Nh /J" dtQ p (t), with C the capacitance of the 
voltage probe. In Refs. [TT] and [H] the FCS of charge fluctuations in a chaotic cavity 
were investigated, conceptually this is similar to the FCS of <j> p addressed here. As 
pointed out in Ref. [11], the quantity J Q dtQ p (t) has no direct physical meaning, but 
can be understood as the time spent by all electrons in the cavity (or probe) after time 
r, or alternatively as the phase picked up due to the potential integrated over time 
r. In analogy to this, the phase n p is proportional to the charge in an energy interval 
dE integrated during r. We only consider the case that there is a large number of 
states in this energy interval. This corresponds to the case of a long delay time tj, of 
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particles inside the probe. For the stochastic path integral approach discussed here 
it is important that the charge on the probe is slowly fluctuating compared to the 
inverse attempt frequency r s = h/eV: 3> t s . This allows to justify the saddle 
point solutions which will be used later on. A different dynamics at voltage probes is 
considered in Ref. [55] . 



3.1. Voltage probe. 

We define a joint probability P(Q, 4> p ) that Q charges are transmitted during r and 
the phase due to the potential of the voltage probe is <p p . By definition we have 
-P(Q) = L d(f> p P(Ci, 4>p). The conjugated variable to <j) p is a number called m p . 
Fourier transforms in the variables Q or 4> p or in both of them are possible. Thus 
there are four equivalent functions containing all informations about the transport 
and the probe statistics: P(Q,0 p ), S(A,m p ), £(A, <j) p ) and ((Q,m p ). Note that the 
counting fields A are periodic in 2-7T, but <p p has a period of 1. Thus the exponent in 
the transformation is multiplied by 2n. Explicitly written we have for example, 

P(Q,<h>) = / > dA^e s ( A ' m ^- 4AQ+2Trl ^ m '', (22) 

nip 

P(Q,M = f dAe^ A ^" lAQ . (23) 



In particular we are interested in the distribution P(4> p ) of the phase <pp alone. It is 
obtained from the relation 

^p) = E p (Q'« = e?(A=o,0p) - ( 24 ) 

Q 

Refs. [2H1 [23 used the stochastic path integral to calculate the generating function 
of a conductor coupled to a voltage or dephasing probe. With this formalism also 
the joint functions P, S, £ and C introduced above are accessible. A description of the 
procedure is given in the appendix, here we simply state the results. 

Most directly the function £(A, <p p ) can be determined by an integral 



C(A,0 p ) = / dXpe s v (A,\ p y p =v<p p )^ ( 25 ) 



A very long measurement time r > defines the stationary case considered here 
where the variables \ p and V p are time independent. The function »Sy(A, X P ,V P ) is 
given by [M] 



Sv(A,X p ,V p ) = T - I dEH (2(») 



with 



Hq = In dct 



1 + n ( \ f S r \S - 1 



(27) 



For a conductor with single mode contacts to the M terminals and to the probe, 
the scattering matrix S is of dimension (M + 1) x (M + 1). The matrix n contains 
the occupation numbers of the different terminals with n = diag(ni, n2, ■ ■ ■ , «m, %>) 
(here n p = f p (V p )), and the matrix A introduces the counting fields, A = 
diag(e iAl , e iX2 , . . . , e iX " , e iX p ) . 

To proceed it is useful to introduce the expressions qm containing multi-particle 
scattering probabilities multiplied with the appropriate counting fields in the contacts. 
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The index I = 0, 1 denotes the number of particles injected into the probe, and k = 0, 1 
the number of particles emitted from the probe [S7J . Below, for the examples of interest 
here, we give the explicite expressions for the q^i ■ With the help of the qj~i , the function 
H is 

H (\, rip) = In [(1 - np) (g o + <7oie iAp ) + n p (q n + gioe - ^)] . (28) 
Using Eq. ||SSJ| we obtain for Eq. $2M, with N = SZV, 

S V =N (<p p ln[q u + q w e- iX »] + (1 - 4> p ) hx[q QQ + <Zoie iA "]) . (29) 
The integral (|25|) is solved in saddle point approximation. The saddle point equation 



^ L = is quadratic in e lXp . The solution satisfying the normalization condition 
/# p e« *) = 1 is 



&i \ v _ -goigio(l - 20p) + y(goigio(l - 20 p )) 2 + 4googoigiogii0 p (l - <j) p ) 

2gio9n(l ^ <A P ) 
The generating function is in saddle point approximation 

e £(A,0 p ) _ e Sv(A,A p ,V0p)^ ^21) 

As demonstrated in Ref . [27J 26] , the saddle point approximation is correct when 
the delay time of the probe is much longer than the inverse average attempt frequency, 
r d > h/eV. 

To obtain the cumulants of the distribution function it is important to look at 
the joint generating function S(A,m p ). The auto- and cross-correlations studied in 
section [2] can be expressed in analogy to Eqs. ([19]) as 

V 2 t d 2 S 
At: 2 dm 2 ' 



=V*(*#) = -?r£™ (32) 



(AV p AI a ) = - (A<fi p AQ a ) - ^- J^ ■ (33) 

The generating function S(A,m p ) is in principle obtained by a Fourier transform 
of £(A, (f> p ) as indicated in Eqs. (|22p and (|2^]) . Another, more transparent way is 
to express 5*(A, m p ) with the help of Eq. ([2"5]) and solve the coupled saddle point 
equations. This is shown in the appendix. 



3.2. Dephasing probe. 

The distribution function for the average occupation number n p of the dephasing probe 
can be found in a way similar to what is described above for the voltage probe. Here, 
the distribution Pe(Q) is energy dependent and the vector Q denotes the charges per 
energy interval, Q = Q(E)dE. The function Pe(Q, n p ) is the joint probability that Q 
charges are transmitted during r and the average occupation number takes the value 
rip in the energy interval dE. The conjugated variable to n p is a number called 7 P , 
and the corresponding generating functions and distribution functions are Pg;(Q,np), 
Se(A,j p ), £e(A, fip) and Ce(Q>7p) w i tn for example 

PB(Q,Tip)= f dAj2 eSE{Anp) ~ lAQ+l2 ™ P ~' P - (34) 

J 7p 
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With this the fluctuations and correlations of section [2] are 

/a ?\ dErd 2 S , , 

^ ="4^' (35) 



(An p Az a ) = t^t^I^- (36) 



e a 2 s 

Note that the counting fields A Q are dimensionless and of order 1, but the variable 7 P 
and the generating function S scale with dEr/h, leading to the proportionality factors 
above. Here again, we are most interested in the distribution function 

P E (n P ) = J2 P E(Q,n P ) = e^ A =°^l (37) 
Q 

The function £e(A, n p ) is given from the stochastic path integral approach in the 
stationary limit: 

e fe(A,n p ) = f dXpe S E (A,\ p ,n p =n p ) = J H (\ p ,n p ) _ (gg) 

The saddle point equation = is quadratic and the solution satisfying the 

normalization is 



e a P= ftoflp (39) 

W g i(l ~ n p ) 



Explicitly we obtain the generating function 



£e(A, rip) = ^^ln 



qoo(l - n p ) + qnnp + 2^J q 10 q 01 n p (l - n p ) 



.(40) 



Even though the distribution function P(Q) of the transmitted charge only is equal 
for both voltage and dephasing probes (in the single channel case), the generating 
functions Eqs. (|3"Tj) and (j4H)l arc different. Technically this originates from the 
different structure of the saddle point equations and their solutions Eqs. and ((501) . 
Physically the differences are based on the different occupation functions of the voltage 
and dephasing probe, a Fermi function and a non-equilibrium occupation function. As 
shown in Fig. [1] the later is a two-step function and carries a probabilistic uncertainty 
leading to the different statistics of the number n p compared to the phase <p p . This is 
already seen in the fluctuations (AV^ 2 ) and (An 2 ) obtained by the Langevin approach 
in section [21 



4. Examples. 

With the formalism presented above, it is possible to discuss the full counting statistics 
of voltage and occupation number fluctuations and their correlations with transferred 
charge of the examples presented in section [2] Here we concentrate on the beam 
splitter, the Mach-Zehnder interferometer and the double barrier. 



4--1- Beamsplitter. 

The beam splitter with transmission T coupled to a probe is shown in Fig. [2] It is 
straightforward to evaluate the equations presented above with the expressions for 
qki ■ The qui are expansion coefficients of the determinant in Eq. (|27p , however for this 
simple example they can be found by considering the different processes in and out of 
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a) b) 




-1.2 -1.2 



Figure 7. Distribution functions P(<p p ) of the integrated voltage (a) and 
Pe{^p) of the occupation number (b) of the the probe of the beam splitter for 
two transmission probabilities. The function P((j> p ) of the integrated voltage 
is a binomial distribution whereas the function P£;(n p ) of the time integrated 
occupation number of the dephasing probe is non-binomial. The plots are 
normalized with respect to N = Si^-L or dN = 



the probe, goo describes the process that no particle enters or leaves the probe. For the 
beam splitter this happens if the particle is reflected at the barrier. This process occurs 
with probability 1 — T. When the particle passes through the barrier with probability 
T, it enters the probe, this is goi- The remaining qiQ and q\\ represent two-particle 
processes: no particle is injected into the probe (probability 1 — T) but one particle is 
emitted from the probe into terminal 2 (counting factor e lX2 ), or one particle enters 
with probability T and one leaves the probe into terminal 2. In conclusion, 



9oo = 1 - T, q 01 = T, 

9io = (l-T)e iA2 , gn = Te^ 2 . 



Inserting this in Eqs. (f3"Tj) and (|4T)| we obtain 

/ Te iX2 1 — T \ 

£(A 2 ,0 P ) =N\^ p ln— + (l-^)ln__J 



£e(A 2 ,™ p ) = diVln 



(l-T)(l-n p ) + JTe^n. 



(41) 

(42) 
(43) 



with N = z¥-L and dN — ^jr- The two functions have a very different form. The 
distribution functions P{(/>p) — e^ ^"^ and PE(n p ) = e^(°'™ p ) are plotted in Fig. [7] 
The function P{4> p ) is a binomial distribution. The saddle point approximation does 
not lead to an binomial factor in Eq. (|4"2"|) . The factor can be obtained by inserting the 
above expressions Eq. (|4Tj) into Eq. (|29|) and evaluating the integral Eq. (|25|) exactly. 
The distribution function of n p for the dephasing probe is however non-binomial as 
can be clearly seen in the figure. 

For this simple example the joint distribution function P(Q,(f> p ) can be directly 
obtained with P(Q,4> P ) = J ^-e^ 2 ^-' 1 ^ . We find that the distribution of the 
phase 4> p equals the distribution of the charge transmitted into a contact of the beam 
splitter: 

P(Q, p ) - eW>M5 Q<h> = P(t p )6 Q4>p . (44) 

Also for the dephasing probe the joint distribution PE(Q,n p ) can be calculated (in 
this function Q denotes the charge per energy interval). The probability function 
PE{Q, n P ) depends on both variables Q and n p , in contrast to the voltage probe, 
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Figure 8. The joint distribution function Pe(Q >n p ) of the charge transferred 
into contact 2 and the integrated occupation number at the dephasing probe of the 
beam splitter for T = 0.5. In contrast for the voltage probe, the joint distribution 
P(Q,4> P ) is a one-dimensional binomial function with P(Q,<f> p ) = P{4>p)&Q,<f> p - 



Eq. (|44p where Q and 4> p act effectively as one variable. This reflects the fact that the 
occupation number h p represents the probability for the emission of a particle, while 
the occupation in the voltage probe determines exactly each scattering event. A plot 
of Pe(Q, n P ) is shown in Fig. [5] 

Although the two generating functions (|4"12|) and (|4"5|) have a very different form, 
integration over <p p and n p respectively leads to the same generating function for charge 
transmitted through a beam splitter: In L d(j) p e^ X2 '^ — J dEln f dn p e^ E ^ X2 ' n ^ — 
N ln[l+T(e* A2 — 1)] . This means that a voltage or a dephasing probe, perfectly coupled 
to the conductor as shown in Fig. [5] does not affect the statistics of transmitted charge. 



4-2. Mach-Zehnder interferometer. 

The Mach-Zehnder interferometer is shown in Fig. [3l The coefficients qui are 
determined by the scattering matrix and are derived in Ref. [27] : 

goo - ((RaRb(1 -e) + T A T B )e lX3 + (R A T B (1 - e) + T A R B )e iXi 

+2^/R A T A R B T B {l - e) cos$(e^ - e lA3 )) (45) 

qai = eR A (46) 
Qio = eT A e l ^+ x ^ (47) 
q n = ({RaRb + T A T B (1 - e))e lX ' + (R A T B + T A R B (1 - e))e lX * 

+2V 'R A T A R B T B {\ - e) cos<$>(e tXi - e^)) (48) 



Using these parameters, the generating functions Eqs. (f31j) and (|40|) become very 
long expressions, which are not explicitely written out here. Plots of the functions at 
A = are shown in Fig. [9j Note that both Eqs. ([31) and (I40|) contain all information 
about the cumulants of the transmitted charge and about the phases <j) p and n p 
respectively. Thus both functions contain terms proportional to cos$ due to the 
interference. However the phase-dependence $ drops out when we study cumulants 




Figure 9. The distribution functions P(<f> p ) (left) of the integrated voltage and 
Pe^p) (right) °f the integrated occupation number of the probe of the Mach- 
Zehnder interferometer for different coupling strength e. The case of symmetric 
beamsplitters is shown in the upper figures and the case of beamsplitters with 
transmission Ta = Tb = 0.7 is shown in the lower figures. For e = 0.1 both 
distribution functions look very similar (in fact they are equal to first order in e). 
For stronger coupling, the tails of the distribution are broader for the dephasing 
probe. The maximum of both distribution functions lies at Ra = 1 — Ta- 



of the phase 4> p or the occupation number n p of the probe. This has been already 
observed in section [2731 and is explained by the setup of the MZI: particles enter the 
probe before they could interfere and no oscillating terms appear in the cumulants of 
the probe. 

Interestingly the distribution functions of <p p and n p given by e^°'^ p ) and e^ B ^ 0,n >'' 
coincide to first order in the coupling parameter s. This means that when only very 
few particles enter the probe, the distribution functions of the potential and of the 
occupation function per energy respectively behave similarly, while for a higher number 
of charges in the probe the differences become more and more clear, in particular the 
function P{<j) p ) is narrower than Pg(n p ). 



4-3. Double barrier. 

The double barrier is shown in Fig. [4] In this setup the electrons enter the probe after 
they interfere and the distribution function depends on the phase picked up. The 
coefficients qu are 



Ri + (1 - e){R 2 + TiT 2 e^) + 2 y /R 1 R 2 (l-e) cos$ d 

<Zoo — ; (49) 

1 + R x R 2 (l - e) + 2 v /i?ii? 2 (l - e) cos$ d 

9oi = = (50) 

1 + i?ii? 2 (l - e) + 2^i? 1 i? 2 (l-e)cos $ d 

9io = Tl£ ; (51) 

1 + R X R 2 {1 - b) + 2^i?ii? 2 (l-e)cos $ d 
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Figure 10. The distribution function P(<j>p) for the integrated voltage (a) and 
Pe (%> ) f° r the integrated occupation number (b) of the probe of the double barrier 
geometry for equal transmission probabilities Ti = T2 = 1/2 and e = 1 — e~ 1 as 
a function of the phase Qj. The different properties of the two probes are most 
apparent around = 7r. 



Plots of the distribution functions of the voltage probe and the dephasing probe, 
Eqs. (|3Tj) and |38|) at A2 = are shown in Fig. [TO] Both functions depend on the 
phase and are most narrow at §d = ft- Around this value the differences between 
the voltage and dephasing probe become visible: the distribution function P((j> p ) is 
narrower than P_E(n p ). For values further away from it the functions look pretty much 
alike. Since the transmission probability T p \ has a maximum in = n, this is in 
agreement with our findings that the more particles enter the probe the less broad are 
the distribution functions. Higher order correlations in the double barrier are briefly 
discussed in the appendix | Appendix B| 

5. Conclusions. 

In this work, using the Langevin approach, we examined the auto-correlations 
of voltage fluctuations at voltage probes and occupation number fluctuations at 
dephasing probes. We investigated correlations between voltage (or occupation 
number) fluctuations at probes and currents at a contact of the conductor. We 
determined these fluctuation spectra for several examples. Subsequently we extended 
this discussion to include higher order cumulants. To achieve this we extended a 
stochastic path integral approach to find the generation function for the distribution 
of voltage and occupation number fluctuations and also the joint distribution of voltage 
(or occupation number) fluctuations and current fluctuations. Interestingly we find 
differences between the fluctuations of voltage and dephasing probes despite the fact 
that for single channel probes the full counting statistics of the transmitted currents 
is identical [551 HI] • This is a consequence of the fact that occupation numbers at 
different energies in the dephasing probe are uncorrelated while in a voltage probe 
they are correlated. The fluctuations of the occupation number of a dephasing probe 
are larger than the potential fluctuations of a voltage probe. For small coupling 
of the probe -corresponding to a small number of particles entering the probe- the 
differences become less important. It is expected that dephasing and which path 
detection are closely related. To understand this we analyzed the sensitivity of the 



qn = 



Ri{l ~e) + R 2 + T x T 2 e iX ^ + 2^/R 1 R 2 (l - e) cos$ d 



(52) 



1 + ^^2(1 - e) + 2 y /R 1 R 2 (l - e) cos $ d 
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correlation functions of voltage (or occupation number fluctuations) and currents to 
Aharonov-Bohm oscillations. For the examples investigated here, we found that the 
Aharonov-Bohm effect in the loop directly probed by the voltage or dephasing contact 
is indeed absent. However, this is a consequence of the chirality of the setups chosen, it 
is not a general feature [53] . Further research is needed to systematically characterize 
the systems where current- voltage correlations are independent of the Aharonov-Bohm 
phase. 
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Appendix A. Stochastic path integral. 

We show how the Eq. (f2"5|) for the voltage probe is obtained starting from the stochastic 
path integral approach [3J [T2] . The voltage probe is described by a Fermi function 
with a potential V p (t) fluctuating on the time scale r^. As derived in detail in Ref. |27J, 
the generating function Sy (A) (the index V stands for voltage probe) is given by the 
stochastic path integral 



ft 



SV(A) = / VV p V\ p e Sv{A ^^ ] (A.l) 



S v = ~ / <H 



-iT d X p eV p + / dEH 



(A.2) 



The function Hq is defined by Eq. (|27|) and is determined by the scattering matrix. 
The term —iTd\ p eV p results from the charge conservation on the probe. The path 
integral in V p and X p is evaluated via a saddle point approximation which is a 
good approximation for t c i 3> h/eV. For a discussion of this important point see 
Refs. [26l[27]. For a long measurement time r 3> T<j, the function Sy is stationary in 
the variables A p and V p , and the saddle point equations are: 

§^ = 0, |^ = ^S v = T -[dEH Q . (A.3) 
d\p dV p h J 

Starting from Eq. (jA.lj) we also find the full counting statistics of the probe. 

We define a joint probability P(Q, <f> v ) that Q charges are transmitted and the 
phase due to the potential at the probe is <f> p . With Eq. (jA.lj) the joint probability 
P(Q, (j) v ) is given by 

P{QA P )= JvV p 5 Up~^ £dtV p (t)j J dAV\ p e^< A ^-y^-^.{kA) 

Here the constraint cj) p ~ 1/(Vt) J dtV p (t) is introduced via a functional delta function 
^ (0p ~ Vt Jo dtv p( 1 )) = Em p ex P faim p ((f> p - f Q T dtV p {t))). Inserting the sum 
we find 

P(Q, <t> p ) = J dA ^ e s ( A > m p)- 4AC !+ 2T ™^ m p (A.5) 
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where the generating function is given by 

e S(A,m p ) = f VVpVXpe s v (A,x p ,v p )-^£-£dtv p (t)_ ( A g) 



In the stationary limit the path integrals in A p and V p are reduced to ordinary integrals. 
Then the sum in m p in Eq. (|A.5|) can be evaluated giving a delta function 5(<fi p — V p /V) 
and thus 

P(Q,4>p) = J dA J dA J , e Sv(A,A P ,V J1 =V* |1 )-iAQ = J dAe «(A,0 p )- l AQ ^ ^ 

leading to Eq. ([25]). 

In order to obtain the joint generating function e s '( A > m p) directly, it is convenient 
to start from Eq. (|A.6|) . In the stationary case, we have to solve the coupled saddle 
point equations in (p p and X p . Up to a constant, the solution is 



S{A,m p ) = iVln 



Qoo + que N P + J (q m - 1n e j + Maiqwe *r 



(A.. 



The cumulants like the auto- and cross-correlations shown in Eqs. (f3"2")l and (|3"3")1 are 
then given by simple derivatives of the above function e s ( A > m p) taken at A = m p = 0. 

For the dephasing probe, the procedure outlined here is exactly the same, but 
all functions concern an energy interval of width dE. However, the joint generating 
function Se(A,j p ) can not be obtained exactly unlike for the voltage probe. 



Appendix B. Higher order correlations. 

In section [5J we found that in all interfering structures considered the current-voltage 
correlations were independent of the phase directly adjacent to the probe. An 
interesting question is whether the full counting statistics discussed in sections [3] and 
[4] gives additional insight on interfering effects in higher order correlations. These 

correlations are joint cumulants like (AV^Ai^) which are proportional to d rfr fi^{i" P ' > 

with k, I € N. With the expression (1A.8|) we can evaluate cumulants of higher order 
for different setups connected to voltage probes. 

For the double barrier setup we find that 5(0, m p ) which determines the 
distribution function of the time-integrated voltage P(<fi P ) depends on the phase in 
the dot 4>d as already seen in section 14.31 All cumulants linear in AV P but to any 
order in the current in lead 2 are independent of the phase <pd and the coupling 
strength e (the proportionality sign hides factors of 27ri): 



<95(A 2 ,ra p ) 



dm p 



e lA2 (l - Rf + AR+ y/e 2i ^{l - i?) 4 + 4e^R(l - R) 2 
00 2(e^(l -R) 2 + 4R) ' ( ' 



The above function at A 2 = determines the mean voltage V p , the first derivative 
gives the current-voltage correlations that were already obtained in section 12.31 etc. 
Cumulants quadratic in AV P however always contain a contribution proportional to 
cos d>d, 



d 2 S(X 2 ,m p ) 



dm 2 



R((e^ (1 - R) + 2R){2 - e) + 4iVT^icos $ d ) 
ro=0 ° C Ne{l-R)e i ^/ 2 {e i ^{l-R) 2 +mf/ 2 ' [ 



For the Mach-Zehnder interferometer, the function 5(0, m p ) depends only on the 
transmission probability of the first beam splitter, Ta, and the coupling parameter 
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e in agreement with our findings of sections 12.31 and 14.21 Also for this example the 
cumulants linear in AV p but to any order in current do not depend on the Aharonov- 
Bohm phase $ or on the coupling strength, but on both T A and T B . We omit to print 
the whole expression which is lengthy. Remarkably, higher order correlations contain 
interfering terms, for example 

d 3 S(A,m p ) ^ R A T A 



dm 2 p d\ 3 Ne 



({R A - T A )(R B - T fl )(2 - e) + A^R A T A R B T B {1 - e) cos$ 



This is an example how full counting statistics can give additional information which 
is not present in the auto- and cross correlations. 
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